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A.  ORIGINAL  OBJECTIVES: 

1 .  Collaborate  with  with  Prof.  Bialynicki-Birula,  who  has  developed  a  "twistorial 
Fourier  transform"  useful  for  generating  EM  and  other  conformal  fields.  We  have 
recently  discovered  a  close  connection  between  this  transform  and  EM  wavelets  and 
intend  using  this  in  several  ways,  including  the  construction  of  twisted  EMW  which 
are  expected  to  have  better  collimation  and  stability  characteristics  than  the  present 
ones. 

2. 1  have  recently  participated  in  a  NATO  Singular  Optics  workshop,  where  I  become 
interested  in  polarization  singularities  and  their  applications  to  EM  wavelets.  I  plan 
to  study  related  issues  of  angular  momentum,  spin,  helicity  and  vorticity.  Some  of 
this  work  may  also  be  done  jointly  with  Prof.  Bialynicki-Birula  as  it  is  related  to  the 
twisted  beams. 

3.  Work  with  Prof.  Konstantin  Lukin,  a  Noise  Radar  expert,  I  met  at  the  above 

workshop,  who  is  interested  in  testing  experimentally  various  anterma  designs  for 
launching  EMW  using  the  sources  I  have  computed. 

4. 1  plan  to  participate  in  the  Multi-scale  Geometry  and  Analysis  in  High  Dimensions 
(MGA)  program  at  IPAM  to  prepare  for  constructing  fast  EM  wavelet  transforms 
based  on  the  recently  computed  Fourier  representation  of  EMW  sources. 

5.  Design  fast  EM  wavelet  transforms  using  ideas  and  methods  from  MGA  such  as 
curvelets,  possibly  with  some  other  participants  from  the  MGA  program. 
OBJECTIVES  FOR  ADDITIONAL  WORK  EFFORT: 

1.  Finish  the  theoretical  analysis  of  PBW  antennas  by  constructing  EM  wavelets  whose 
polarization  conforms  to  the  geometry  defined  by  the  propagation  of  the  oblate 
spheroidal  wave  fronts.  The  current  formulation  is  based  on  Whittaker  potentials. 


which  have  no  relation  to  the  spheroidal  geometry.  Finding  conformal  polarizations  has 
practical  consequences  since  it  would  greatly  simplify  fhe  mafhemafical  expressions  for 
fhe  sources,  hence  fheir  implemenfafion,  and  should  also  improve  fhe  efficiency  of  fhe 
anfennas. 

2.  Collaborafe  wifh  Professors  Anfhony  Devaney  and  Edwin  Marengo  (Norfheasfern 
Universify)  fo  help  Dr.  Richard  Albanese  and  his  feam  af  Brooks  Cify  Base  consfrucf 
realizafions  of  PBW  wifh  optimal  power,  energy  and  efficiency  for  fhe  purpose  of  high- 
resolufion  remofe  sensing  and  imaging. 

B.  OBJECTIVE  ADDED  DURING  GRANT  PERIOD: 

Eind  nafural  polarizations  for  EM  wavelefs  fhaf  are  consisfenf  wifh  fhe  spheroidal 
geomefry  of  fhe  scalar  pulsed-beam  wavelefs  (PBW)  used  fo  derive  fhe  EM  wavelefs. 

STATUS  OE  EEEORT 

[A  brief  statement  of  progress  towards  achieving  the  research  objectives.  Please  make  this  substantive 
(Limit  to  200  words).] 

1 .  Prof.  Bialinicki-Birula's  mefhod  fumed  ouf  fo  be  capable  only  of  giving  sourceless  EM 
waves,  which  rules  ouf  ifs  applicabilify  fo  fhe  production  of  EM  wavelefs.  However,  I 
developed  an  independenf  mefhod  for  consfrucfing  EM  wavelefs  wifh  arbifrary 
angular  momenfum  abouf  fhe  beam  axis.  See  fhe  affached  papers  Retarded  multipole 
wavelets  Parfs  I  and  II. 

2.  Prof.  Lukin  visifed  me  in  July  2005  and  his  efforfs  culminafed  in  a  reporf  (affached), 
buf  Dr.  Albanese  and  I  concluded  fhis  had  limifed  value. 

3.  The  collaboration  wifh  Profs.  Devaney  and  Marengo  resulfed  in  a  significanf  paper 
(affached),  fo  appear  in  IEEE  Transactions  in  Anfennas  &  Propagafion,  which 
validafes  fhe  realizabilify  and  usefulness  of  EM  wavelefs. 

4.  The  mefhod  of  Debye  pofenfials  (suggesfed  by  Prof.  Devaney)  yields  EM  wavelefs 
whose  polarizations  are  consisfenf  wifh  fhe  spheroidal  geomefry  of  fhe  scalar  PBW.  I 
believe  all  fhe  elemenfs  are  now  in  place  for  fhe  realizafion  and  fabrication  of  EM 
wavelef  sources. 

ACCOMPLISHMENTS /NEW  EINDINGS 

[Describe  research  highlights,  their  signiticance  to  the  field,  their  relationship  to  the  original  goals,  their 
relevance  to  the  AF's  mission,  and  their  potential  applications  to  AF  and  civilian  technology  challenges.] 

1 .  The  singular  pulsed-beam  wavelet  sources,  supported  on  disks,  have  been  replaced 
by  equivalent  shell  sources  supported  on  spheroidal  surfaces  or  shells  surrounding  the 
disks. 


2.  This  method  was  recently  found  to  yield  a  conserved  transmission  current  which 
expresses  the  invariance  of  fhe  fransmission  amplifude  befween  wavelef  sources  wifh 
respecf  fo  deformation  of  fhe  equivalenf  source. 
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Beginning  in  October  2007, 1  have  been  giving  3-hour  monthly  lectures  on  EM  wavelets 
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Abstract 

Electromagnetic  wavelets  are  constructed  using  scalar  wavelets  as 
superpotentials,  together  with  an  appropriate  polarization.  It  is  shown  that 
oblate  spheroidal  antennas,  which  are  ideal  for  their  production  and  reception, 
can  be  made  by  deforming  and  merging  two  branch  cuts.  This  determines  a 
unique  field  on  the  interior  of  the  spheroid  which  gives  the  boundary  conditions 
for  the  surface  charge-current  density  necessary  to  radiate  the  wavelets.  These 
sources  are  computed,  including  the  impulse  response  of  the  antenna. 

PACS  numbers:  02.30.Jr,  02.30.Uu,  41.20.Jb,  41.85.Ct 


1.  Complex  distance  and  its  branch  cuts 

Electromagnetic  wavelets  were  introduced  in  [K94]  as  localized  solutions  of  Maxwell’s 
equations.  They  are  ‘wavelets’  in  the  historical  sense  of  Huygens  as  well  as  the  modern 
one:  being  generated  from  a  single  ‘mother  wavelet’  by  conformal  transformations  including 
translations  and  scaling,  they  form  frames  that  provide  analysis-synthesis  schemes  for  general 
electromagnetic  waves.  Together  with  their  scalar  (acoustic)  counterparts,  they  have  been 
called  physical  wavelets.  It  was  pointed  out  that  they  can,  in  principle,  be  emitted  and  absorbed 
causally,  and  applications  to  radar  and  communications  have  been  proposed  [K96,  K97,  Kl] 
based  on  their  remarkable  ability  to  focus  sharply  and  without  sidelobes.  Similar  objects, 
long  studied  in  the  engineering  literature  under  the  name  complex-source  pulsed  beams,  have 
been  used  to  build  beam  summation  methods  and  analyse  the  behaviour  of  solutions  (see 
[HEl]  for  a  recent  review).  However,  to  implement  the  proposed  applications  to  radar  and 
communications,  the  wavelets  must  be  realized  by  simulating  their  sources.  This  has  proved 
to  be  difficult,  and  detailed  investigations  have  only  been  made  recently  [HLKO,  K3].  Here  I 
present  a  new  and  rather  complete  analysis  of  the  sources  based  on  an  insight  which  I  believe 
is  a  key  to  their  realization:  the  sources  can  be  constructed  from  the  very  branch  cuts  that  give 
the  wavelets  their  remarkable  properties. 


0305-4470/04/225929H- 1 9$30.00  ©  2004  lOP  Publishing  Ltd  Printed  in  the  UK 


5929 


Physical  wavelets  are  based  on  the  idea  of  displacing  a  point  source  to  complex 
coordinates.  Since  a  real  translation  gives  nothing  new,  it  suffices  to  discuss  a  point  source 
with  purely  imaginary  coordinates  ia.  It  will  be  seen  that  this  results  in  a  real,  coherent, 
extended  source  distribution  parametrized  by  the  single  vector  a,  much  as  an  antenna  dish  can 
be  described  by  a  single  vector  giving  the  orientation  and  radius  of  the  dish. 

Recall  the  definition  of  the  complex  distance  a  from  the  imaginary  source  point  ia  to  the 
real  held  point  r, 

a(r  —  ia)  =  \/  (r  —  ia)  ■  (r  —  ia)  =  \l r'^  —  —  2ia  •  r.  (1) 

For  each  hxed  source  location  ia  ^  0,  its  branch  points  form  a  circle: 

CT  =  0=j>reC  =  {re]R^  :  r  —  a,  a-r  —  0}. 

It  is  important  to  note  the  topological  difference  between  a  ^  0,  when  —  C  is  multiply 
connected,  and  a  =  0,  when  C  contracts  to  the  origin  and  —  C  becomes  simply  connected. 
Writing 

a  ^  p-iq  (2) 

note  that  ( 1 )  implies 

2  2  2  2 

r  —  a  —  p  —  q  a  -  r  —  pq 

from  which  one  easily  obtains  the  relations  to  the  cylindrical  coordinates  with  z-axis  parallel 
to  a: 

az  —  pq  ap  —  a\/ r'^  —  z}  —  si  P^  +  a^\l  a^  —  q^.  (3) 

This  gives  an  important  bound  on  the  ‘complexness’  of  a  in  terms  of  the  ‘complexness’  of  its 
argument: 


\q\  <  a 

It  follows  from  (3)  that 


or  |ImCT(2:)|  ^  |Im2;| 


z  —  r  —  la. 


2  2  2  2 
P  Z  _  P  Z_ 

p2  q2  _  qi  q2 

hence  the  level  surfaces  of  p^  form  a  family  of  oblate  spheroids 
Sp  —  {r  :  p^  —  const  >  0}  = 


r  : 


2  2 
P  Z 


-  +  —  =  1 
+  p^  p^ 


and  those  of  q^  form  the  orthogonal  family  of  one-sheeted  hyperboloids 


Tiq  =  {r  :  0  <  q^  —  const  <  a^j  = 


a^-q 


2  q^ 


To  complete  the  picture,  we  include  the  following  degenerate  members  of  these  families. 


(4) 

(5) 


5o  =  {r  :  0  ^  p  ^  fl,  z  =  0}  (6) 

Hq  —  {r  a  ^  p  <  oo,  z  —  Q] 

Ha  —  {r  ■.  p  —  Q,  — oo  <  z  <  oo)  (7) 


where  Sq  and  Hq  each  form  a  twofold  cover  of  the  indicated  sets. 

The  families  Sp  and  Tiq  are  depicted  in  figure  1,  along  with  the  azimuthal  half-plane 
<p  —  constant.  They  are  confocal,  with  the  branch  circle  C  as  their  common  focal  set.  Note 
that  the  intersection  of  Sp  with  Hq  consists  of  two  circles  whose  further  intersection  with  the 


Figure  1.  The  level  surfaces  of  a  foiTn  an  oblate  spheroidal  coordinate  system. 


azimuthal  half-plane  consist  of  two  points  for  each  choice  of  ,  0).  When  p  —  q  —  Q, 

the  two  circles  merge  with  the  branch  circle  C.  The  set  of  numbers 

(a,  <p)  =  (p,  q,(p)  —oo<p<oo  —a^q^^a  0  ^  (p  <  Ijt 

therefore  gives  a  twofold  cover  of  —  C.  To  obtain  a  coordinate  system,  we  must  choose 
between  the  two  covers,  and  this  amounts  to  choosing  a  branch  cut  that  makes  a  single  valued, 
as  explained  below.  This  will  result  in  a  one-to-one  correspondence  between  (ct,  </>)  and  points 
r  e  not  on  the  branch  cut,  giving  an  oblate  spheroidal  coordinate  system. 

If  we  continue  cr  analytically  around  a  closed  loop  threading  the  branch  circle  C,  it  returns 
with  its  sign  reversed.  To  make  it  single  valued,  it  is  therefore  necessary  to  choose  a  branch 
cut  that  prevents  the  completion  of  the  loop.  Note  from  (1)  that 

r  ^  a  ^  a  ^  ±{r  —  iacosO).  (8) 

The  spatial  region  with  r  a  will  be  called  the.  far  zone  (we  need  a  >  0  here  to  set  the  scale). 
Since  we  want  cr  to  generalize  the  usual  positive  distance  r,  we  insist  that 

r^a^p^r  q^acosO.  (9) 

If  follows  from  (9)  that  the  branch  cut  B  is  bounded  since  it  must  be  entirely  contained  inside 
any  spheroid  Sp  with  p^  sufficiently  large,  and  its  boundary  must  be  the  branch  circle: 

3B  =  C.  (10) 

(The  alternative  is  a  branch  cut  extending  from  C  to  inhnity,  but  this  violates  (9).)  B  is  therefore 
a  membrane  spanning  C,  and  any  such  membrane  will  do.  The  situation  is  best  understood 
topologically.  The  analytic  continuation  of  the  distance  has  opened  up  a  window  connecting 
the  two  branches  of  r  =  ^/r  •  r,  thus  making  —  C  multiply  connected.  The  spherical 
coordinates  r  and  9  merge  analytically  into  a ,  which  is  double  valued,  and  the  choice  of  a 
branch  cut  B  makes  B  simply  connected  and  a  single  valued. 

Let  CTo  denote  the  complex  distance  with  the  flat  disc  So  as  branch  cut: 

cro  =  Po  -  igo  Po^O  -aifqoi^a.  (11) 

The  complex  distance  ag  with  B  as  branch  cut  is  now  defined  as  follows.  Choose  an  arbitrary 
reference  point  tq  in  the  far  zone,  where  ag  =  ctq-  To  find  erg  at  any  other  point  r,  continue 
analytically  along  an  arbitrary  path  from  tq  to  r,  with  the  following  rule:  whenever  the  path 
crosses  B,  ag  changes  sign.  This  gives  a  unique  definition  of  erg,  and  both  pg  and  qg  have  a 
jump  discontinuity  across  the  interior  of  B.  Of  course,  erg  =  ctq  =  0  on  C. 


Figure  2.  The  cut  function  /  {q)  =  Xe  {q)  with  s  =  0.005  and  its  branch  cut  with  a  =  \.  The  two 
sheets  have  been  purposely  separated  to  show  the  double  cover. 

A  general  branch  cut  can  be  specified  by  a  function  /  (p,  (p)  as 

B  —  {r  \  p  —  X{q,<l>),  —a  ^  q  ^  a,0  ^  (j)  ^  2tx],  (12) 

where  the  cut  function  x  must  satisfy 

x(-q,(p)  =  -x(q^4>)  xiq,^^)  =  x(q,0).  (13) 

The  first  condition  ensures  that  a  changes  sign  across  B,  while  the  second  ensures  that  it 
is  continuous  across  the  half-plane  0  =  0  on  each  side  of  B.  Note  that  B  need  not  be 
cylindrically  symmetric.  We  will  be  especially  interested  in  the  cuts  Ba  defined  by  the 
cylindrical  cut  functions 

X„(g)  =  Q;Sgn(g).  (14) 

Note  that  on  Be,  we  have 

az  —  pq  —  aq  Sgn(g)  —  a\q\. 

If  a  >  0,  then  the  values  q  0  generate  the  upper  spheroid 
{r  eSa  :  z>  0}. 

But  this  does  not  include  the  branch  circle  C,  so  the  bounding  condition  (10)  is  not  satisfied. 
The  problem  is  that  Sgn(^)  is  undefined  at  q  =  0,  and  the  set  of  all  points  with  ^  =  0  is 
the  degenerate  hyperboloid  Ho  (7).  Hence  we  define  the  part  of  Bo,  with  q  —  0  as  the  apron 
bridging  the  gap  between  C  and 

=  {r  :  fl  <  p  ^  \/a2~+~a2,  z  =  0). 

Thus,  for  a  >  0,  Xaiq)  defines  the  upper  spheroidal  branch  cut 

=  (15) 

Similarly,  x-aiq)  defines  the  lower  spheroidal  cut 

dB.e,=C.  (16) 

As  O'  ^  0,  both  cuts  contract  to  the  doubly  covered  flat  disc  spanning  C,  which  is  the 
degenerate  spheroid  So  (6). 

Every  branch  cut  B  is  doubly  covered.  Consider  any  simply  connected,  closed  surface  iS 
containing  C  in  its  interior.  Think  of  N  as  a  balloon  and  of  C  as  a  rigid  wire  ring.  Now  deflate 
the  balloon,  and  you  have  a  branch  cut  bounded  by  C.  In  particular,  if  we  take  S  —  Sa  and 
keep  the  upper  spheroid  rigid  while  deflating,  the  balloon  stretches  around  the  ring  to  cover 
the  underside  of  5,^  and  we  obtain  the  cut  (15)  (see  figure  2). 

The  image  of  a  branch  cut  as  a  balloon  enclosing  the  singular  ring  is  similar  to  Penrose’s 
idea  of  cosmic  censorship  in  general  relativity  [W99,  chapter  5]  where  a  horizon  (N)  prevents 
the  outside  observer  from  seeing  a  naked  singularity  (C).  In  light  of  the  connection  with 


Newman’s  analytic  Coulomb  field  (see  the  discussion  below  (54)),  the  two  may  in  fact  be 
closely  related. 

The  discontinuity  of  Xa  in  (14)  causes  two  problems;  the  q  —  Q  contribution  is  undefined 
(hence  the  aprons  had  to  be  chosen  ‘by  hand’),  and  the  resulting  cut  had  a  sharp  edge.  For 
computational  purposes,  it  may  be  better  to  use  smooth  cut  functions  to  avoid  both  problems. 
Let  e  >  0  and  define 

=  -Imln^^'*'’;^)  xiq)  ^oiX.iq).  (17) 

JT  \e  —  iq  J 

For  s  a,  Xi.{q)  is  &  smoothed  version  of  Sgn(q')  and  the  resulting  branch  cut  closely 
approximates  Bo,  without  the  need  to  define  the  apron  separately.  This  is  shown  in 
figure  2. 


2.  Scalar  wavelets 


For  any  fixed  choice  of  branch  cut  B,  we  now  denote  the  complex  distance  simply  by  ct.  Scalar 
wavelets  are  then  defined  as  the  retarded  solutions 

—  —  z  —  r  —  ia  x  —  t  —  ib  (18) 

a  a 

where  we  have  set  the  propagation  speed  c  =  1  (otherwise  gr  —  —  cr/c))  and  g  is  the 

analytic-signal  transform  of  a  driving  signal  go(0>  defined'  as  the  convolution  of  go  with  the 
Cauchy  kernel: 


g(T^)  = 


_l_  p  goit')dt' 

2i7T  y_oo  r  -t' 

/•“  goit')dt' 

'2-yt  i_oo  (C  - 1)^  +  b^ 
gi{t,b)-t-ig2(t,b). 


X  —  t  —  ib 

J_  /•“  {f -t)go{t')dt' 

27r  i_oo  (C  -  0^  + 


(19) 


g\(t,  b)  and  g2{t,  b)  are  smoothed  versions  of  ^goit)  and  its  Hilbert  transform,  with  b  as  the 
smoothing  parameter.  We  assume  that  go(0  decays  at  infinity,  from  which  it  follows  that  gix) 
is  analytic  in  the  upper  and  lower  complex  time  half-planes  C*.  The  original  driving  signal 
can  be  recovered  as  the  boundary  value 


go(t)  =  lim  [g{t  -  ib)  -  g{t  +  ib)] . 
i^+O 

(The  limit  of  the  sum  gives  the  Flilbert  transform.)  In  particular,  if  go  vanishes  on  any  open 
interval  /,  this  interval  becomes  a  window  between  the  upper  and  lower  half-planes  through 
which  the  functions  g(f  ±  ib)  can  be  connected  so  that  they  are  both  part  of  a  single  analytic 
function.  (This  is  a  special  case  of  the  edge  of  the  wedge  theorem  in  higher  dimensions;  see 
[K3].)  Since  every  practical  driving  signal  vanishes  at  least  in  the  remote  past,  this  property 
will  be  assumed.  Note  that  this  excludes  time-harmonic  driving  signals,  which  are  however 
idealizations. 

Now  consider  the  numerator  of  (18), 


g{x  -  a)  ^  git  -  p  -  iib  -  q)). 

Suppose  that  \b\  ^  a.  Then  gix  —  a)  is  undefined  along  the  semi-hyperboloid  where 
qir)  —  b,  except  when  t  —  pir)  is  in  the  zero-set  of  go-  On  the  other  hand,  if  \b\  >  a,  then 
b  —  qir)  vanishes  nowhere  and  gix  —  a)  is  analytic  at  all  (r,  t).  Therefore  we  assume  from 
now  on  that 


\b\  >  a 


(20) 


This  is  a  special  case  of  a  multidimensional  definition;  see  [K3]. 


1 


so  that  g(T  —  cr)  is  defined  unambiguously  everywhere.  The  imaginary  source  coordinates 
must  therefore  belong  either  to  iht  future  cone  or  to  the  past  cone  of  spacetime, 

b  >  a  ^  (a,  b)  e  V+  b  <  —a  ^  (a,  b)  e  y_ 


which  means  that  the  complex  4-vector  from  the  source  point  iy  =  i(a,  h)  to  the  field  point 
X  —  {r,  t)  belongs  either  to  the  forward  tube  or  the  backward  tube  of  complex  spacetime 
[SW64,  K3], 

T±  —  {z  —  x  —  iye'C^\ye  V±}  x  —  {r,t)  y  =  (a,  b).  (21) 


The  source  distribution  of  'T  is  now  defined  as  a  generalized  function  by  applying  the  wave 
operator, 

5(z)  =  (9,2  -  v2)4/(z)  =  a 4/(2),  (22) 


where  indicates  that  the  operator  acts  only  on  the  real  spacetime  variables  x  of  the  field 
point.  It  is  well  known  that 

h(t  —  r) 

□  — - ^  =  4jTh(t)8(r)  (23) 

r 

for  any  differentiable  function  h,  and  this  can  be  extended  to  'h(z).  Since  'T  is  differentiable 
in  r  everywhere  outside  of  the  branch  cut  B,  (23)  suggests  that  5  is  a  (Schwartz)  distribution 
supported  on  B,  a  conclusion  borne  out  by  a  rigorous  analysis  [K3].  The  discontinuity  of  4' 
across  B  gives  a  combination  of  simple  and  double  layer  terms  of  S  on  B  [K3]. 

The  frequency  content  of  g{x)  determines  that  of  4^  and  should  therefore  be  understood. 
Substituting  the  Fourier  representation  of  go  into  the  definition  (19)  and  reversing  the  order  of 
integration  give 


1  /■“  df'  r 

=  T-  /  - / 

2i7r  r  -t'  J_, 


dco^  If 

i_oo  2:T^°^‘^^2i7r  J_ 


2n 

oo 


(24) 


The  contour  in  the  second  integral  can  be  closed  in  the  lower  half-plane  if  h  >  0  and  in  the 
upper  half-plane  if  h  <  0,  giving 

/OO 

— go(«)0(«&)  r  =  r  -  ih  (25) 

■OO 


where  0  (mb)  is  the  Heaviside  step  function.  Thus  if  b  >  0,  g  contains  only  the  positive- 
frequency  components  of  go,  and  if  h  <  0,  it  contains  only  the  negative-frequency  components. 
In  either  case,  the  factor  e“"*  in  the  extended  Fourier  kernel 

ia>T 

e  =  e  e 


acts  as  a  low-pass  filter,  substantially  damping  frequencies  \oj\  ;:§>  and  thus  smoothing  out 
g(x).  If  the  driving  signal  go  is  assumed  real,  then  go(±<w)  are  related  by  complex  conjugation 
and  therefore  so  are  g(f  =F  tb).  If  go  is  complex,  then  go(i<^)  are  unrelated  and  so  are 
g(t  =F  i^)- 


Example.  Let  g(x)  be  the  (n  —  l)st  derivative  of  the  Cauchy  kernel^, 


g(T)  =  C„(r)  =  (i9,)"-';^ 

ziTrr 


(n  -  1)! 

27ri"r" 


(26) 


^  The  driving  signal  is  the  singular  distribution  go(0  =  (i3t)"  *15(1),  but  this  can  be  approximated. 


whose  Fourier  transform  is 


Cnico,  b) 


At  e'“'C„(?  -  ib)  =  Sgn(h)©(®/7)tt."“‘ 


e-"* 


(27) 


Thus,  while  b  acts  to  suppress  high  frequencies,  n  >  1  acts  to  suppress  low  frequencies  and 
we  end  up  with  a  band-pass  filter  whose  effective  centre  frequency  and  bandwidth  are  given 
by  a  Poisson  distribution, 

n  Jn 

«„  =  —  Ao)  =  - .  (28) 

b  \b\ 

The  behaviour  of  T*  in  the  far  zone  is  governed  by  that  of  g(r  —  s).  By  (9), 

C  (T  -  a)  =  («-!)! _ 

"  2:/ri"(T— ct)"  27t[(b  —  a  cosd) -t- i(t  —  r)]" 

is  a  pulse  with  angle-dependent  duration 


T (9)  —  \b  —  a  cos0|  ^  \b\  —  a  —  T^nm  >  0 


(29) 


being  shortest  at  0  =  0  if  h  >  a  and  at  0  =  tt  if  h  <  —a.  While  the  pulse  duration  is 
independent  of  n,  the  strength  of  the  peak  depends  jointly  on  the  size  of  n  and  the  smallness 
of  b  —  a\ 

MW  =  IS(r (30) 

To  get  a  measure  of  the  diffraction  angle,  assume  h  >  a  for  definiteness.  Fix  fi  >  0  and  look 
for  the  angle  0^  at  which 

MiOft)  =  e"^M(0). 

Then 


(b  —  a  cos0^)"  =  e^(h  —  a)" 

which  gives 

2sin2(0^/2)  =  1  -  cos0^  =  (e^/"  -  (31) 

Thus,  6^  can  be  made  small  either  by  choosing  b  —  a  a,  or  n  ;:§>  fi.  In  either  case,  the 
right-hand  side  gives  0^/2.  A  reasonable  measure  is  obtained  with  /S  =  1 . 


3.  From  scalar  to  vector  wavelets 


It  is  well  known  that  every  electromagnetic  field  can  be  derived  from  a  pair  of  real 
scalar  potentials,  the  most  well-known  examples  of  which  are  the  Whittaker  and  Debye 
superpotentials  [N55].  In  this  section  we  use  the  scalar  wavelet  'T  as  a  complex  Whittaker 
superpotential.  Although  this  is  equivalent  to  using  a  pair  of  real  potentials,  disentangling 
the  real  and  imaginary  parts  leads  to  unnecessarily  complicated  expressions,  something  like 
taking  the  real  and  imaginary  parts  of  a  complicated  analytic  function  iy)  in  order  to 

obtain  two  real  harmonic  functions.  To  see  how  bad  it  gets,  note  from  (49)  that  the  fields  and 
currents  contain  terms  of  the  type  —  a)  with  k  =  0,  1,  2.  In  the  simplest  case  k  —  2 

(which  will  give  the  radiation  terms  of  the  field),  (19)  gives 

gjT:  -  g)  _  gi{t  -  b,b  -  q)  +  ig2{t  -b,b-q) 
a  P  ~  iq 


Re 


g(^  -  g) 
a 


pgijt  -b,b-q)-  qg2it  -b,b-q) 
p^  +  q^ 


(32) 


Im 


g(^  -  g) 
a 


Pgiit  -b,b-  q)-i-qgiit  -b,b-q) 
p^  +  q^ 


and  it  is  clear  that  the  real  expressions  quickly  become  unmanageable.  Thus,  although  we 
work  with  complex  potentials  and  fields,  we  view  this  as  a  very  compact  and  efficient  way 
of  computing  the  real  fields.  In  particular,  our  expressions  contain  nothing  extraneous  since 
their  imaginary  as  well  as  real  parts  have  a  direct  physical  significance.  This  strategy  is  based 
on  the  analyticity  of 'T  outside  the  source  region,  which  will  indeed  make  harmonic  pairs  out 
the  fields  D  and  B,  as  seen  below. 

With  'T  as  a  complex  Whittaker  superpotential,  we  define  the  retarded  complex  Hertz 
potential 

Z  =  T'Tr  (33) 

where  vr  is  a  fixed  complex  polarization  vector  that  can  be  seen  [K3]  to  be  a  combination  of 
(real)  electric  and  magnetic  dipole  moments.  The  real  and  imaginary  parts  of  Z 

Z  —  Ze+  iZ,„ 

are  interpreted  as  electric  and  magnetic  Hertz  vectors  [BW99,  p  84,  85].  They  generate  a 
4- vector  potential  (/x  =  0,  1,  2,  3)  by 

Ao  =  -V-Z,  A  =  9,Z,  +  VxZ,„  (34) 

which  automatically  satisfies  the  Lorenz  condition 

3,Ao  +  V-A  =  0.  (35) 

In  turn,  it  follows  from  potential  theory  (or  the  Poincare  lemma  for  differential  forms)  that 
every  4-vector  potential  satisfying  (35)  can  be  written  in  the  form  (34),  so  this  representation 
is  quite  general.  (We  can  even  dispense  with  the  Lorenz  condition  by  performing  a  gauge 
transformation  on  A^.  See  [N55]  for  an  excellent  and  thorough  account  of  Hertz  potentials 
and  their  enormous  gauge  group.)  The  real  vector  fields  Pg  and  defined  by 

P  =  P«-HiP,„  =  DZ  =  (□4')7r  =  57r  (36) 

are  the  electric  and  magnetic  polarization  densities.  They  are  distributions  supported  spatially 
on  the  branch  cut  B.  Since  we  are  in  Lorenz  gauge,  the  charge-current  density  is  7^  =  □  A|U, 
hence 

-Pe  J  =  d,Pe  +  V  X  P„,  (37) 

with  charge  conservation  guaranteed  by  the  Lorenz  condition.  The  polarization  densities  thus 
act  as  ‘potentials’  for  the  charge-current  density,  a  property  inherited  directly  from  (34). 

The  reason  why  Hertz  potentials  will  be  so  useful  can  be  seen  by  computing  the  fields: 

P  =  V  X  A  =  V  X  V  X  Z„,  -I-  9,V  X  Z,  (38) 

and 

E  =  -VAo  -  d,A  =  ■  Ze-  dfZe  -  9,V  x  Z,„ 

= -DZ^-h  V  X  V  X  Z,  -  9,V  X  Z„,.  (39) 

Taking  into  account  (36)  gives 

D^E  +  Pe^V  xV  X  Ze-dtV  X  Z,„  (40) 

which  is  a  kind  of  ‘harmonic  conjugate’  of  (38),  so  the  real  fields  D,  B  can  be  expressed 
compactly  in  the  complex  form  [S41,  pp  32-34] 

P  =  P-i-iP  =  V  X  V  X  ZH-i9,Z.  (41) 

Note  that  outside  the  branch  cut  B,  Pe  —  0  and  D  —  E.  The  Hertz  formalism  thus 
automatically  takes  account  of  the  polarization,  so  that  the  expression  (40),  if  interpreted 
as  a  distribution,  is  valid  even  within  a  singular  source  region. 


Again  I  emphasize  that  F  =  D  +  iB  is  ‘real’  in  the  sense  that  D  and  B  are  real,  physical 
fields.  Yet  F,  like  'h,  is  analytic  in  the  source-free  complex  spacetime  region 

Tb  =  {(r  -ia,t-  ib)  :\b\>a,r^  B}.  (42) 

More  simply,  because  of  their  spheroidal  symmetry,  'k  and  F  are  analytic  functions  of  the  two 
complex  variables  (a,  t)  in  the  region 

Ub  =  {(cr,  t)  e  :  \b\  >  a,  p  ^  Xsiq)]  (43) 

where  xeiq)  is  the  cut  function  for  B.  Thus  D  and  B  are  really  harmonic  conjugates 
as  suggested  earlier.  On  the  other  hand,  P  and  characterize  the  singularities  spoiling 
analyticity  in  the  source  region,  including  the  branch  points  and  branch  cuts.  This  differs 
from  the  usual  practice  in  the  frequency  domain,  where  (D,  B)  are  the  real  parts  of  separate 
complex  fields  (Dc,  Be),  and  it  might  appear  that  these  two  representations  are  in  conflict 
since  the  real  fields  cannot  be  extracted  by  taking  the  real  and  imaginary  parts  of  Dc  +  iBc. 
To  clarify  this,  consider  the  frequency  components  of  F, 

/OO 

■OO 

and  note  that  since  F  is  complex,  its  positive-  and  negative-frequency  components  are 
independent .  Therefore  a  general  monochromatic  field  consists  of  two  terms, 

G™ono  =  e-‘'“'F„  +  e‘“'F_^  «>0 

=  e-‘'“'(D„  +  iB„)  +  e-“'(D_„  +  iB.,) 

=  e-^'d)^  +  iB„)  +  +  iBl) 

where  the  reality  conditions  have  been  used  on  D^,  and  Ba,.  The  representation  of  a 
monochromatic  field  is  therefore  no  different  from  that  of  a  general  held; 

dniono  —  ^mono  tBjjjQjjQ 

with  both  helds  real: 

£>mono(0  =  2Re{e-‘“'£>„}  B^.UO  =  2Re{e--“'B„}. 

Recall  that  for  b  >  a  and  b  <  —a,  the  analytic  signal  g(r)  contains  only  positive-  and 
negative-frequency  components.  Therefore 

b  >  a  ^  Deo  +  tBio  =  0Vft)<0  b  <  —a  Du,  +  iB®  =  0  V  w  >  0. 
This  shows  that  the  monochromatic  components  of  the  electromagnetic  wavelets  satisfy 

Bmono(0  =  2Re{ie-‘“"r>„}  =  D  mono  TtlleS), 

SO  B  trails  D  if  b  >  a,  and  it  leads  D  if  b  <  —a.  (Recall  also  that  the  pulse  travels  along 
±a  if  ±b  >  a.)  More  generally,  the  wavelets  are  helicity  eigenstates  with  helicity  1  if  h  >  a 
and  —  1  if  h  <  —a.  This  concept  applies  not  only  the  time-harmonic  components  but  also 
to  general  time  domain  helds  [K3a].  As  already  mentioned,  using  the  analytic  combinations 
of  helds  also  has  the  great  advantage  of  compactness  and  simplicity  over  the  alternative  of 
disentangling  the  real  and  imaginary  parts. 

Before  launching  into  the  held  computations,  I  want  to  prepare  the  way  for  computing 
equivalent  currents  on  the  spheroid  Sa  in  the  coming  section.  Let  us  construct  this  spheroid 
from  the  two  branch  cuts  B±a  given  in  (15)  and  (16).  Let  us  now  denote  by  a  the  complex 
distance  with  the  disc  Sq  as  branch  cut.  This  will  be  used  as  a  reference  for  dehning  the 
complex  distance  functions  with  B±a  as  cuts,  which  we  denote  by  a±.  Let  V±  be  the  volumes 
bounded  by  together  with  So,  so  that 

3y+  =  5^-50  9y_  =  <So-5- 


<?■_  =  <7 


Figure  3.  Values  of  the  branches  (T±  of  the  complex  distance  function  determined  by  the  branch 
cuts  B±a,  given  in  terms  of  the  branch  a  determined  by  the  disc  Sq. 


where  the  signs  are  related  to  the  orientations  of  iS*  and  Sq  hy  a.  The  union  and  compliments 
of  V±  will  he  denoted  by 

V  =  v+uv-  vX  =  -  y±  y'  =  -  y. 


Now  recall  the  rule  for  crossing  a  branch  cut  other  than  the  reference  cut  5o:  erg  changes  sign. 
Thus,  denoting  the  complex  distance  functions  with  respect  to  B±a  by  ct±,  we  have 


a  in  y^ 

—a  in  V± 


(44) 


as  shown  in  hgure  3. 

The  held  radiated  /oi/if/y  by  the  two  branch  cuts  B±a  is  therefore 


'2F(a,T) 

in 

r 

F„(cr,  t)  = 

F{—a,  r)  +  F{a,  r) 

in 

F{a,  r)  +  F(—a,  r) 

in 

y_. 

Observe  that  there  is  no  field  discontinuity  in  going  from  V+ 

to  y_ 

2F(a,  t) 

in 

y' 

Fa(a,  x)  = 

F{—a,  r)  +  F{a,  r) 

in 

y 

(45) 


as  depicted  in  hgure  4. 

The  transition  a  —>■  —a  across  a  branch  cut  turns  retarded  helds  into  advanced  helds 
since 


'T(— cr,  r)  —  — 


.g(7  +  cr) 


(46) 


a 

Although  advanced  helds  are  usually  associated  with  acausal  behaviour,  there  is  a  perfectly 
causal  explanation  for  (46).  Consider  the  held  radiated  backward  from  B^,  as  observed  in 
y+.  Due  to  the  curvature  of  the  back  side  of  Be  this  held  converges  towards  the  focal  ring  C 
and,  having  passed  though,  it  is  no  longer  in  V+  and  therefore  diverges  normally.  A  similar 
argument  explains  why  the  held  emitied  forward  from  hrst  converges  towards  C  and 
then  diverges  away  from  it.  The  usual  acausal  behaviour  associated  with  advanced  helds  is 
due  to  the  assumption  that  they  remain  advanced  for  the  indehnite  future.  (This  argument 
also  applies  to  time-reversed  acoustics  [FO],  where  time  reversal  occurs  only  in  a  bounded 
spacetime  region.) 


Figure  4.  Interior  and  exterior  fields  radiated  by  the  oblate  spheroid  Sa,  represented  as  a 
combination  of  the  two  branch  cuts  B±a  • 


It  was  shown  in  [K3]  that  the  sources  of  t)  and  a,  r)  are  equal  and  opposite; 
that  is  they  form  a  source-sink  pair: 

□  'hl-a,  r)  = r)  = -5.  (47) 

The  proof  is  trivial  for  real  point  sources,  where 

□  ^  ±  r)&{r)  =  -Ango{t)&{r). 

r 

But  it  is  more  subtle  for  complex  point  sources  because  the  extended  delta  function 
,  1 

l(z)  — -  z  —  r  —  \a 

Ana 

with  a  7^  0  fixed,  is  not  supported  at  a  single  point  but  on  the  entire  branch  cut  B  and  therefore 
f{a)Uz)  ^  fiO)S(z). 

In  fact,  the  left-hand  side  is  not  even  defined  since  cr  is  discontinuous  precisely  on  the  disc 
supporting  S{z)',  therefore,  some  care  must  be  used  in  proving  (47). 

Equation  (47)  shows  that  the  interior  superpotential  'T(a,  r)  +  'T(— cr,  r)  is  sourceless, 
as  are  the  Hertz  potentials  and  electromagnetic  fields  derived  from  it.  The  interior  field  is 

Fo(a,  r)  —  F(a,  r)  +  F(—a,  r).  (48) 

Let  us  first  compute  the  exterior  field  F,  which  will  give  the  interior  field  by  symmetrizing 
with  respect  to  a .  Let 

I'' ^  d„‘l‘ — ^  g,  -  S,g,(T  -  a) 

a  cr^  a 

z 

u  =  Va  =  V/7  —  iVq  =  — 
a 

and  note  that  it  is  a  complex  unit  vector: 

2  •  Z 

U  •U  =  - T-  =  1. 

<7^ 


Thus 


V  X  Z  =  VT*  X  TT  =  'T'm  X  TT 

V  X  V  X  Z  =  X  (m  X  tt)  +  'T'V  x  (it  x  tt) 

and  by  a  simple  computation, 

M  X  (it  X  Tt)  =  Xu  —  TV  X  —  U-TV 

Xu  +  TV 

V  X  {u  X  tv)  — - . 

a 

Therefore  (41)  gives 

'T' 

F  —  ^>"{Xu  —  tv) - {Xu  +  Tv)  +  i'T'ii  X  vr 

a 


F= I 


1  i  Sr  Sr  Sr 
Xu-\  —  +  —r  +  — 
a 


a  a 

This  expression  will  be  written  compactly  as 
F  —  LXu  —  Mtv  —  iNu  X  TV 

where 


t,^. 


3gr 


r  _  Sr  _^^Sr  '^Sr 

Lr  +  «  +  ,2 

G  G^  G^ 


Sr 


a 


u  X  TV.  (49) 


(50) 


+  (51) 

a 


We  now  examine  the  far  field  to  see  under  what  conditions  the  polarization  vector  tv  gives  the 
strongest  beams.  In  the  far  zone  (9)  we  have 

r^a^a'^r  —  ia  cos  6  u  ^  By. 

Therefore 


-F’far  = 


g(t  -  cr) 


(XBr  —  TV  —  IBy  X  Tv) 


g(T  -  cr) 


(tTx  +  IBy  X  TTx) 


(52) 


where 


TV±_  —  TV  —  {tv  •  er)ey 

is  the  component  of  tv  orthogonal  to  r  which,  as  expected,  is  the  only  one  that  matters  in  the 
far  zone.  Note  that  while  we  have  replaced  cr  by  r  in  the  denominator  of  (52),  the  presence 
of  Im  a  —r  cos  0  in  g(r  —  cr)  plays  an  essential  role  in  determining  both  the  collimation  of 
the  beam  and  the  duration  of  the  pulse,  as  already  seen  in  (30)  for  g{x)  —  C„{r). 

The  far  field  satisfies  the  helicity  condition 

IBy  X  Ffc  =  Ffai-  (53) 

or  equivalently 

J5far  =  X  JDfar  Df^  =  -By  X  Bf^y. 

As  we  are  interested  mainly  in  the  paraxial  region  of  the  far  zone,  the  most  efficient  choice  of 
TV  is  orthogonal  to  a.  Since 

r  —  ia  p  z  —  ia 


r  —  p  +  za  ^  u 


la  p  z 

a  a  a 


this  implies 


X  —  U  •  TV  — 


p  •  TV 


The  far-zone  energy  density  is 


far  I 


•  +  |Sf„|2}  =  i|Ffc|2 


and,  by  (53),  the  far-zone  Poynting  vector  is 


*^far  —  -^far  ^  -f^far  —  -^far  ^  -^far  —  t 

2l 

=  X  (er  X  Ffar)  =  ^|Ffar|^e, 


far  X  -Ffar 
—  ^far^r 


since  Ffai  is  orthogonal  to  . 


4.  Equivalent  currents 


In  principle,  the  scalar  source  generates  the  charge-current  density  by  (36)  and  (37).  But  this 
would  involve  not  only  the  messy  disentangling  of  the  real  and  imaginary  parts  of  P  =  Stv 
(with  both  factors  complex),  but  also  dealing  with  the  singular  nature  of  S.  While  S  is  well 
defined  mathematically  as  a  distribution  [K3],  it  seems  to  be  of  little  direct  value  from  a  practical 
point  of  view.  Since  S  is  supported  on  the  branch  cut  B,  one  expects  the  electromagnetic  sources 
to  consist  of  a  surface  charge  density  jo  and  a  surface  current  density  j.  But  it  turns  out  that 
these  surface  sources  are  singular  on  the  branch  circle  C,  where  ct  =  0.  The  essence  of  the 
problem  can  be  understood  from  a  careful  analysis,  given  in  [Kla],  of  a  much  simpler  case, 
which  we  now  recall. 


Example.  The  analytically  continued  Coulomb  field  due  to  point  charge  of  strength  g  =  1  is 
„  1  r  —  ia 

C(r-ia)  =  -V- - =  - - (54) 

47tG 


Newman  [N73]  has  shown  that  this  can  be  identified  with  a  real  electromagnetic  field  (Z),  B) 
by 

C^D  +  iB  (55) 


interpreted  as  the  flat-spacetime  (zero-mass)  limit  of  the  Maxwell  field  in  the  Kerr-Newman 
solution  in  general  relativity  [N65],  which  represents  a  spinning  black  hole  of  unit  charge^. 
It  is  instructive  to  compute  the  surface  sources  on  a  branch  cut,  which  for  simplicity  we  now 
take  to  be  the  disc  So  defined  in  (6).  On  the  upper  and  lower  faces  of  5o,  we  have 

p  ^  +0  z  ^  ±0  a  — >•  =pi-y/a2  _  p2 


hence 


C  ^  =F- 


ip 


47T{a^  —  p2)^/2  ^  47T(a^  — 

The  jumps  in  D  and  B  across  the  cut  are  therefore 

a  _  2p 


&D  = 


27T(a^  —  p2)3/2 


SB  = 


2Tt{a^  — 


Since  SD  is  orthogonal  and  SB  is  tangent  to  So,  it  follows  that  the  magnetic  suiface 
charge-  and  current-densities  vanish  as  required.  The  electric  surface  densities  are  given  by 
[J99,p  18] 


jo  =  —a  ■  SD  —  — 
j  —  a  X  SB  —  — 


2jT(a^  —  p^)^/2 
ca  X  p 


2jt{a^  —  p^)^/2  2jt{a'^  —  p^)^/^ 


(56) 


3  When  D  in  (55)  is  reinterpreted  as  aNewtonian  force  field,  then  B  is  a  gravitomagnetic  field  related  to  the  ‘dragging’ 
of  Einsteinian  spacetime  in  the  vicinity  of  a  spinning  body.  Evidently,  this  effect  survives  the  flat-spacetime  limit  as 
the  conjugate-harmonic  partner  to  the  Newtonian  gravitation. 


where  we  have  inserted  the  speed  of  light  (taken  earlier  to  he  c  =  1)  for  dimensional  reasons. 
Before  discussing  the  problem  with  (56),  note  that  if  we  define  the  local  charge  velocity  by 


3  cp 

V  =  —  =  —60 

a  a 


(57) 


its  linearity  in  p  suggests  a  ‘hydrodynamic’  interpretation  of  Sq  as  a  rigidly  spinning  charged 
disc  with  angular  velocity 


=  c/a.  (58) 

In  particular,  the  rim  C  is  moving  at  the  speed  of  light.  While  this  conclusion  seems  bizarre  in 
ordinary  electrodynamics,  it  is  entirely  consistent  with  the  origin  of  the  field  C  as  the  residual 
Maxwell  field  of  a  charged,  spinning  black  hole.  The  investigation  in  [Kla]  has  sparked 
a  renewed  interest  in  Newman’s  original  paper  [N73],  leading  to  similar  interpretations 
of  linearized  gravitational  helds  [N2]  and  a  generalized  Lienard-Wiechert  field  where  the 
radiating  point  source  moves  along  an  arbitrary  trajectory  in  complex  spacetime  [N4].  Our 
antennas  will  be  similar,  but  their  source  is  a  dipole  following  a  complex  trajectory  and  not  a 
monopole,  and  so  their  charge-current  densities  are  generated  by  polarizations. 

We  now  come  to  the  main  lesson  taught  by  this  example.  C  is  an  analytic  continuation 
of  the  Coulomb  field  of  a  point  source  with  charge  <2=1-  If  the  continuation  is  to  make 
physical  sense,  the  total  charge  should  remain  unchanged.  This  is  contradicted  by  jo,  which 
is  not  only  strictly  negative  but  whose  total  charge  on  iSo  is  —  oo!  To  resolve  this  difficulty, 
it  is  necessary  to  treat  the  charge-current  density  as  a  singular  volume  distribution,  just  as  the 
scalar  source  S  =  □  T'  was  treated  in  [K3].  The  inhomogeneous  Maxwell  equations  now 
state  that  the  (volume)  charge-  and  current-density  are 

To  =  V  •  F  J  =  -drF  -  iV  X  F  (59) 


while  the  homogeneous  Maxwell  equations  require  that  Jq  and  J  be  real.  Taken  as  definitions 
of  the  sources  in  the  sense  of  generalized  functions,  it  was  shown  in  [K3a]  that  (59)  indeed 
give  a  sensible  answer.  The  equivalent  surface  sources  on  any  spheroid  Sa  with  a  >  0  are 
defined  by 

Jo  =  ep-SF  j  =  -iBp  X  SF  (60) 

where  the  outgoing  unit  normal  on  Sa  is  computed  in  the  appendix.  These  sources  are 
found  to  be  complex,  which  means  that  they  include  a  magnetic  charge-current  density;  the 
latter  vanishes  in  the  limit  a  — >•  0,  in  agreement  with  the  above  conclusion.  The  advantage 
of  using  a  >  0  is  that  the  sources  are  smooth  and  hounded,  with  a  total  charge  Q  =  1  as 
required.  As  a  0,  they  decompose  into  surface  sources  on  the  interior  of  the  disc  which 
coincide  with  (56),  plus  line  sources  on  the  rim  C.  The  line  sources  carry  a  total  charge  of  oo, 
but  when  the  entire  source  distribution  is  treated  as  a  generalized  function,  it  carries  the  correct 
total  charge  Q  =  1.  The  problem  with  (56)  is  that  the  jump  conditions  (using  infinitesimal 
pillboxes  and  loops)  can  be  applied  only  on  the  interior  of  the  disc  and  not  on  its  boundary 
C.  A  similar  argument  applies  to  every  branch  cut,  showing  that  caution  must  be  exercised  in 
computing  equivalent  sources,  a  lesson  we  will  recall  when  computing  the  currents  required 
to  produce  electromagnetic  wavelets. 

Finally,  we  turn  to  computing  the  equivalent  sources  for  F  on  the  spheroid  Sa-  Some 
important  properties  of  equivalent  real  scalar  surface  sources  were  analysed  in  [HLKO],  but 
their  connection  to  the  vector  case  and,  specifically,  to  our  topological  use  of  branch  cuts, 
remains  to  be  explored. 

According  to  (45)  and  (48),  the  jump  in  the  field  across  the  spheroid  is 
&F{a,  t)  =  F(cr,  t)  -  F(-cr,  t) 


(61) 


where  the  complex  distance  a  with  respect  to  iSq  is  continuous  across  Sa-  Unlike  the  sum  (48) 
of  retarded  and  advanced  fields,  the  difference  (61)  does  have  sources  and  they  are  confined  to 
the  surface  Sa,  which  we  shall  presently  compute.  Begin  by  writing  (50)  in  the  more  explicit 
form 


F{a,  r)  =  L{a,  r)ku  —  M{a,  x)-k  —  iN{a,  x)u  x  tt 
with  L,  M,  N  given  in  terms  of  the  retarded  signal  gr{o,  x)  —  g{x  —  a)  by 


(62) 


Tt  ^  Sr  ^^gr  ^  3gr 
L{g,  t)  =  —  +  ^  ^ 

G  G^  G^ 


A/fr  \  8r  8r  8r 

Mia,  x)  =  —  +  —r  +  — 

G  G^  G^ 


AT/  \  8r  8r 
Nia,  x)^  —  +  —r. 

G  (T^ 


Define  the  mixed  signals  g±  by 

g±{a,  r)  =  gix  -a)±  gix  +  a) 


and  note  that 


a  —a  ^  u  — - —u  X  —  u  •  TV  ^  —X. 

a 

Then  we  obtain  the  following  expression  for  the  field  discontinuity: 

&F  —  L{a,  x)Xu  —  Mia,  x)tv  —  iNia,  x)u  x  tv 


(63) 


where 

Lia,  t)  =  —  +  +  — 

G  G^  G^ 


Mia,  r)  =  ^  ^  ^ 

G  G^  G^ 


AT/  \ 

Nia,  t)  =  —  +  — . 

a  a^ 


Before  going  on  to  compute  the  currents,  note  that  (45)  can  be  modified  so  that  the  interior 
field  is  any  source-free  solution  of  Maxwell’s  equations,  i.e., 

IFia,  r)  in  V' 

Fint(r,r)  in  U 

V  ■  Fint  =  0  ifirFint  =  V  X  Fint. 


F  = 

±  Ci  - 


The  choice  of  an  interior  solution  other  than  Fo(cr,  r)  will,  of  course,  modify  the  equivalent 
sources  on  Sa-  However,  unless  Fint  fits  into  the  spheroidal  geometry,  the  resulting  sources 
can  be  expected  to  be  much  more  complicated  and  unnatural,  and  probably  will  not  benefit 
from  the  ‘magic’  of  complex  source  points.  Probably  the  most  general  class  of  interior  fields 
that  do  fit  the  geometry  consists  of  arbitrary  multiples  of  Fq,  i.e. 


Fintir,  t)  =  vFoia,  x). 


Then  (61)  is  replaced  by 

SFia,  x)  —  pFia,  x)  —  uF(— a,  r)  pL  +  v  =  2  (64) 

and  all  computations  below  easily  generalize  to  this  case.  However,  only  when  pv  —  v  —  \ 
can  the  radiating  surface  he  interpreted  as  a  combination  of  branch  cutsl  I  believe  that  this 
case  is  the  most  natural  and  expect  it  also  to  be  the  most  useful.  For  this  reason,  only  it  will 
be  treated  here,  although  our  results  easily  extend  to  the  case  (64). 

I  now  compute  or  estimate  the  various  inner  and  outer  products  needed  in  (60).  Free  use 
will  be  made  of  the  results  derived  in  the  appendix,  and  there  is  no  pretense  of  rigour.  I  will 
assume  that 


0  <  a  <$;  a 

which  means  that  the  spheroid  is  rather^at.  By  (3)  and  (75), 

y - 

p'^  f  a^  —  Z  —  aq 


|Vp|  — 
O' 


Z  —  aq 


Thus 


r  —  ia  p+  {z  —  ia)a  p  —  ia 

u  —  - =  - - .  (65) 

a  a  a 

Recall  that  n  is  orthogonal  to  a,  so  that 


TT  —  TtpEp  +  Jt^e^ 


and  thus 

r  •  TT  p-TV  pjXp 

X  —  U  •  TV  —  -  =  -  =  - 

a  a  a 

From  (76)  in  the  appendix,  the  outgoing  unit  normal  on  Sa  is 


Vp  ar  +  qa  ^  qa  q  ^ 

e„  = - =  .  ^  ,  =  — a.  (66) 

\'^  p\  ^  a?-  +  q^-J  +  q'^  Icrl 

The  term  has  been  retained  in  the  denominator  to  control  the  singularity  at  the  equator. 
(This  is  the  main  advantage  of  using  Sa  instead  of  Sq.)  The  approximation  (66)  fails  very  near 
the  equator  q  —  Q,  where  is  far  from  parallel  to  a,  but  the  analysis  in  [HLKO]  suggests  that, 
for  scalar  wavelets  at  least,  the  immediate  vicinity  of  ^  =  0  can  be  ignored.  More  precisely, 
it  was  shown  that  for  time-harmonic  driving  signals  of  frequency  co,  the  effective  aperture, 
emitting  most  of  the  radiation,  consists  of  the  front  surface  of  the  disc  So  parametrized  by 

k~^  ^  q  ^  a  i.e.  p^^a^—l/k^  k  —  eojc.  {61) 


Of  course,  this  has  significance  only  if  >  1 .  Lower  frequencies  generate  mostly  a  reactive 
field  that  swirls  around  the  source  region  and  is  eventually  reabsorbed^.  Thus,  to  obtain  a 
high  radiation  efficiency,  it  is  necessary  to  use  signals  g(t)  with  little  low-frequency  content, 
such  as  linear  combinations  of  high-order  derivatives  of  the  Cauchy  kernel  (26).  (Of  course,  a 
careful  repetition  of  the  analysis  needs  to  be  made  specifically  for  the  electromagnetic  case.) 

The  inner  products  needed  to  find  Jo  are 


ep-u  =  ep-  (Vp  -  iVq)  =  |Vp|  Rs  — 

pI 

q  „ 

Ep  ■  TV  ^  - a  ■  TV  —  0 

I 

Ep  ■  {u  y.  tv)  —  TV  ■  (Ep  y  u)  —  TV  ■  (Ep  y  (V p  —  iVq)) 

ip 

=  i|Vq'|7r  ■  (Eg  X  Ep)  =  i|V^|7r  ■  60  ]— r^0- 


The  outer  products  needed  for  j  are 


ep  X  16  =  i|Vg|e0  Rs 

q  „  iff 

e„  X  TT  — a  y  TV  ^  — (rtnE^  —  jrj,p) 
p  IcrL  ^  ^ 


Ep  y  (u  y  tv)  —  (Ep  ■  Tv)u  —  (Ep  ■  u)tv 


-{Ep  ■  U)tV  RS - TV. 


Using  these  in  (60)  gives  the  approximate  surface  charge  density 


a\(j\jQ  R»  Lapjtp  ■¥  Napjv^ 


(68) 


^  This  is  consistent  with  the  general  fact  that  ‘DC  components  do  not  propagate’.  It  is  also  the  basis  of  one  of 
the  close  connections  between  electromagnetic  wavelets  and  mathematical  wavelet  theory,  since  it  amounts  to  an 
admissibility  condition  on  electromagnetic  wavelets  [K94,  p  214]. 


and  the  approximate  surface  current  density 


o'lcrlj  ^  (Lp^TTp  —  Ma^iCp  +  NapTt^)e^  +  +  Naanp)ep.  (69) 


As  expected  from  our  example  of  the  analytic  Coulomb  potential,  the  equivalent  sources  on 
a  spheroid  are  complex,  indicating  the  presence  of  unrealizable  magnetic  charges.  Since  the 
magnetic  sources  in  that  example  vanished  as  a  — >•  0,  it  is  reasonable  to  hope  that  this  will 
also  be  the  case  here.  As  the  spheroid  Sa  with  0  <  a  a  is  very  flat,  it  may  be  possible 
to  choose  the  phase  of  the  polarization  vector  tt  (representing  the  mixture  of  electric  and 
magnetic  dipoles)  so  as  to  minimize  the  magnetic  sources  over  Sa  excluding  the  vicinity  of 
the  rim  C,  and  the  latter  region  can  be  ignored  for  highly  oscillatory  driving  signals  as  shown 
in  [HLKO].  This  question  will  be  addressed  in  detail  elsewhere. 

Finally,  we  compute  the  impulse  response  of  the  antenna,  i.e.  the  sources  when  the 
driving  signal  is  the  impulse 

goit)  =  Sit)  ^  sir)  =  — ^  =  Ci(t). 

ZIJTT 


Note  that  the  real  point  source  version  of  the  scalar  wavelet  (1 8)  is  then  the  retarded  propagator 
for  the  wave  equation. 

Sit  —  r) 

'h  ^  'I'o('‘,  t)  —  -  □  fl'o  =  Alt  Sir,  t)  (70) 

r 

where  the  precise  relation  between  and  fl'o  is  given  in  terms  of  complex-distance  potential 
theory  in  [K3].  The  mixed  signals  are 

1  It  7  7 

g+  — - 1 - =  -  where  u  —  r  —  a 

2i7r(T  —  ct)  2i7r(r-i-CT)  mu 

1  1  cr 

^  2i7r(r  —  ct)  2i:7r(T-i-CT)  inu 


and  their  time  derivatives  are 

r^  +  CT^ 

8+  ;  ^ 
2a  r 


8-  = 


mu^ 


This  gives 

IScr^r  -  IOct^t^ -I- 3r5 

L  =  - - 


8+  = 


8-  = 


2t^  +  6tct^ 
ijTu^ 
6r^CT  -I-  2ct^ 
mu^ 


M  = 


9a‘^r  -2a^r^  +  r 


N  = 


3ct4-i-6ct2t2  -  t^ 


which  can  be  substituted  into  (68)  and  (69)  to  obtain  the  impulse  response. 

In  view  of  the  discussion  following  (67),  we  are  actually  more  interested  in  the  system’s 
response  to  the  handhass  signal  in  (26), 


C„(t)  =  (i3,)”-iCi(T) 


(n  -  1)! 
27ri"r'' 


(-a,)"-'Ci(T). 


The  induced  surface  source  can  be  computed  directly  from  the  impulse  response: 


4"’  =  i-3b) 


tJ- 

n—\ 


JfJ.‘ 


(71) 


5.  Concluding  note 


Source-free  relativistic  fields  always  extend  analytically  to  the  double  tube  domain  7±  (21)  of 
complex  spacetime,  as  explained  in  [K3].  I  find  it  quite  remarkable  that  the  extension  'h  (a,  r) 
of  the  propagator  (70)  generates  fields  with  spatially  compact  sources  that  are  analytic  in  the 
source-free  parts  Tg  of  complex  spacetime  obtained  by  removing  the  world  tubes  swept  out 
by  the  sources.  The  boundary  values  of  these  analytic  fields  then  characterize  the  singular 
sources,  as  shown  above. 
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Appendix 


The  complex  unit  vector  u  is  given  by 


hence 


Note  that 


u  —  Vct  =  —  —  y  p  —  iS/q 
a 


Vp 


pr  +  qa 
p^  +  q^ 


Vq 


pa  —  qr 
p^  +  q^ 


and 


which  gives 


2  =  1 

Vp- 

Vg  =  0 

+  a^ 

9 

P~~ 

q^  +  2a^ 

p^  +  q^ 


p  +  2  +  q^ 


7  9 

2  p  +  a 


-  2^  2 
p^  +  q^ 


2 

p^  +  q^ 


The  unit  vectors  in  the  directions  of  increasing  p  and  q  are  therefore 
pr  -I-  qa  pa  —  qr 


^p^  +  q^^p'^  +  a^ 


y/p^+q^y/a^  -  q^ 


(72) 

(73) 

(74) 


(75) 


(76) 
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Abstract 

In  the  companion  paper,  a  charge-current  distribution  was  obtained  for  radiating 
electromagnetic  wavelets.  However,  it  cannot  be  realized  because  of  two 
drawbacks:  (a)  it  requires  the  existence  of  magnetic  charges,  and  (b)  the  charge- 
current  distribution,  which  is  concentrated  on  a  spheroidal  surface  Sc ,  is  too 
singular  for  practical  implementation.  Both  of  these  difficulties  are  resolved 
here.  The  first  is  resolved  by  using  Hertz  vectors  to  generate  a  charge-current 
distribution  on  Sa  due  solely  to  bound  electric  charges.  The  second  is  resolved 
by  replacing  Sa  with  a  spheroidal  shell  of  finite  thickness.  This  generalizes  the 
usual  boundary  conditions  on  an  interface  between  electromagnetic  media  by 
allowing  the  transition  to  be  gradual. 

PACS  numbers:  02.30.Jr,  02.30.Uu,  41.20.Jb,  41.85.Ct 

(Some  figures  in  this  article  are  in  colour  only  in  the  electronic  version) 


1.  Introduction 

This  paper  is  aimed  at  physicists  and  mathematical  physicists  interested  in  classical 
electrodynamics  and  optics.  In  the  companion  paper  [K4],  a  coherent  charge-current 
distribution  for  radiating  electromagnetic  wavelets  was  constructed  on  an  oblate  spheroidal 
surface  Sa-  Its  main  drawback  was  the  necessity  of  including  magnetic  along  with  electric 
charges,  making  the  sources  impossible  to  realize.  Here  we  show  how  this  difficulty  can  be 
overcome  by  using  Hertz  vectors  [BW99,  J99]  to  generate  a  charge-current  distribution  on 
Sa  due  solely  to  bound  electric  charges.  However,  this  distribution  is  still  too  singular  to 
realize.  We  show  how  it  can  be  replaced  by  a  simple  volume  distribution  on  a  spheroidal 
shell  containing  Sa-  Our  method  generalizes  the  usual  boundary  conditions  on  an  interface 
between  electromagnetic  media  by  allowing  the  transition  to  occur  gradually,  without  incurring 
additional  complexity. 
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We  rely  on  the  concepts  in  [K4],  with  some  improvements  in  the  notation.  For  further 
background  on  physical  wavelets  and  complex-source  pulsed  beams,  see  [K94,  HLKO,  HFl], 
The  complex  distance  from  the  imaginary  source  point  ia  to  the  real  observation  point  r  will 
be  denoted  by 

r  =  (r  —  ia)  •  (r  —  ia)  —  p  —  iq  (1) 

and  the  complex  time  by 

t  —  t  —  ib.  (2) 

The  imaginary  time  h  plays  the  role  of  an  overall  scale  parameter,  similar  to  the  scale  of 
ordinary  wavelets  in  one  dimension,  determining  the  duration  of  the  pulsed-beam  wavelets. 
The  imaginary  space  vector  a  similarly  controls  the  spatial  extent  and  orientation  of  the 
wavelets.  The  real  and  imaginary  parts  of  f  satisfy  the  inequalities 

\p\  ^  r,  \q\  ^  a  where  r  —  |r|,  a  —  |a|.  (3) 

For  fixed  a  7^  0,  the  branch  points  of  f  in  form  the  circle 
C  —  {r  f  —  Q]  —  {r  :  a  •  r  —  0,  r  —  a] 
and  the  ‘standard’  branch  of  f ,  defined  by  p  ^  0,  has  for  its  branch  cut  the  disc 
V  —  {r  \  p  =  0}  =  {r  :  a  •  r  =  0,  r  ^  fl}. 

Here  f  is  real-analytic  in  except  for  a  jump  discontinuity  due  to  a  sign  reversal  upon 
crossing  V.  Every  other  branch  satisfying  the  positivity  condition 

o— 

can  be  obtained  by  continuously  deforming  P  to  a  membrane  B  with  the  same  boundary, 
dB  =  dV  =  C. 


The  associated  branch  of  f  is  defined  by 


fB  =  PB-  iqB 


r, 

r, 


r  ^Vb 

r  e  Vb 


(4) 


where  Vb  is  the  compact  volume  swept  out  by  deforming  V  to  B.  It  follows  [K4a]  that  fg  is 
real-analytic  in  except  for  a  sign  reversal  across  B. 

The  scalar  wavelet'  with  branch  cut  B  is  defined  by 


'kg  — 


gji  -  re) 
rB 


where  g  is  the  analytic  signal’  associated  with  a  driving  signal  g{t)  exciting  the  source  by 


gQ)  = 


j_  /•“  g(t')dt' 
Ini  y_oo  t  —  t' 


(5) 


which  is  indeed  analytic  in  the  complement  of  the  support  of  g: 
dgO)  1 

— ^  = -(9,  -  i9i)g(t  -  ih)  =  0  Vr^suppgcM.  (6) 

at*  2 

The  significance  of  the  extension  parameter  h  can  be  understood  by  noting  that  the  real  and 
imaginary  parts  of  g  are 


gbit)  = 


b  /■“  git')  dt' 

277  it  -  ’ 


gbit) 


J_  /■“  it' -t)git')dt' 
277  it  -  t')^  +  b^ 


*  Strictly  speaking,  the  term  ‘wavelet’  should  be  reserved  for  certain  choices  of  g,  as  explained  in  [K4].  The  scalar 
wavelet  of  order  n  is  obtained  with  g  as  the  /?th  derivative  of  the  Cauchy  kernel  C(i)  =  l/lTTii.  Also,  note  that  we 
use  units  in  which  the  propagation  speed  is  c  =  1 . 


Figure  1.  {p,q,4>)  fonn  an  oblate  spheroidal  coordinate  system. 


Here  gt  is  a  smoothed  version  of  g  with  h  as  the  scale  or  resolution  parameter,  while  gi,  is  a 
smoothed  version  of  the  Hilbert  transform  of  g,  again  with  h  as  the  scale  parameter.  Thus, 
time  variations  of  order  less  than  \b\  are  suppressed  in  g{t  —  ib). 

Due  to  its  denominator,  'hg  is  singular  on  C,  where  fg  =  0,  and  discontinuous  in  the 
interior  of  B,  where  fg  reverses  sign.  To  avoid  any  further  singularities,  we  want  to  ensure 
that  the  numerator  g{t  —  fg)  is  analytic  in  all  of  and  for  this  it  suffices  to  have  its  argument 
bounded  away  from  the  real  axis  by  a  positive  distance.  Since 

f  -  fg  =  t  -  pg  -  i(h  -  qs)  (8) 

and  |(?g|  =  1^1  ^  fl,  a  necessary  and  sufficient  condition  is 

fl  <  \b\.  (9) 


This  states  that  the  imaginary  spacetime  four-vector  (a,  b)  is  time-like,  belonging  to  the  future 
cone  of  spacetime  if  a  <  b  and  the  past  cone  if  h  <  —a.  The  condition  (9)  will  be  assumed 
from  now  on,  making  'kg  real-analytic  in  —  B. 

The  source  Eg  of  'kg  is  defined  by  applying  the  wave  operator: 

47rEg  =  n'kg  where  A.  (10) 

Eg  can  be  easily  shown  to  vanish  wherever  'kg  is  twice  differentiable,  hence 


r  ^  S  ^  Eg  =  0. 


(11) 


To  characterize  the  source  on  B,  we  must  apply  □  in  a  distributional  sense  [GS64].  Just  as 
differentiating  the  Heaviside  function  gives  the  delta-function,  differentiating  a  discontinuous 
function  like  'kg  in  a  distributional  sense  gives  a  single  layer  on  the  surface  of  discontinuity, 
represented  by  a  delta-function  of  a  variable  normal  to  that  surface.  Since  Eg  is  obtained  by 
differentiating  'kg  twice,  it  will  consist  of  a  combination  of  single  and  double  layers  on  B 
[K4,  K4a].  Moreover,  these  layers  diverge  on  the  boundary  C  —  dB  since  'kg  is  singular  there. 
This  singularity  will  be  tamed  below  by  combining  wavelets  with  different  branch  cuts. 

The  variables  (p,q)  defined  by  (1),  together  with  the  azimuthal  angle  0  about  the 
a-axis,  determine  an  oblate  spheroidal  coordinate  system  where  the  level  surfaces  of  p  are  the 
spheroids 


5,: 


x^  +  y^ 

p2  +  qZ 


p^O 


(12) 


and  the  level  surfaces  of  q  are  the  orthogonal  hyperboloids  Hq .  All  these  quadrics  are  confocal, 
having  the  circle  C  as  their  common  focal  set.  This  is  depicted  in  figure  1.  As  p  0,  Sp 
shrinks  to  a  double  cover  of  the  disc  P. 


Figure  2.  The  upper  hemispheroidal  branch  cut  B*  with  its  apron  Aa  ■ 


When  the  source  of  'hg  is  computed,  it  will  be  singular  on  C  due  to  the  singularity  of 
'hg  there.  In  the  case  of  a  real  point  source,  this  corresponds  to  the  singularity  of  S(r)  at 
r  =  0.  Recall  that  the  latter  can  be  regularized  by  replacing  the  origin  with  a  sphere  of  small 
radius  r  —  a,  whence  S(r)  is  replaced  by  a  uniform  distribution  on  the  spherical  surface.  The 
delta-function  can  then  be  dehned  in  terms  of  the  limit  a  — >•  0.  The  equivalent  procedure  now 
is  to  replace  the  sphere  by  the  oblate  spheroid  Sc  which  is  defined  by  p  =  Re  f  =  a  >  0.  But 
we  can  go  a  step  further  and  represent  Sa  as  a  sum  of  two  branch  cuts,  something  that  cannot 
be  done  for  a  real  point  source  since  the  deformation  of  a  point  is  still  a  point.  Thus,  consider 
the  branch  cut 


(13) 


consisting  of  union  of  the  upper  hemispheroid 
5+  =  {r  e  :  z  >  0} 

and  the  apron 

=  {r  :  r  •  a  =  0,  ^  ^  H-  a^} 

connecting  5^  to  the  branch  circle  C,  as  shown  in  figure  2.  (The  apron  must  be  included  so 
that  dB^  =  C  as  required.) 

Similarly,  let 


(14) 


be  the  union  of  the  lower  hemispheroid 
=  {r  e  :  z  <  0} 

with  Ac-  For  simplicity,  denote  the  complex  distance  with  branch  cut  B^  by  r±  instead  of 
rg±  and  the  corresponding  wavelet  by  'I'i.  Let  be  the  interiors  of  the  upper  and  lower 
hemispheroids  and  14  be  the  interior  of  Sa-  According  to  (4), 


r  e  Va  ^  t)  =  'I'(-f,  f),  'I'_(r,  t)  =  'h(r,  t) 
r  e  Va  ^  'h+(r,  f)  =  4'(r,  ?),  'h_(r,  t)  =  '!'(— f,  t) 

r  ^Va  ^  'h+(r,  t)  =  t)  =  4'(r,  t). 


Consider  the  average  of  'I'± , 


'hA(T,  t)  =  5{'I'_(r,  t)  +  4'+(r,  i)}. 


(15) 


Then  by  the  above, 

'TA(r,  t)  = 


t), 


reVc 

riV, 


where  the  internal  field  'I'l  and  the  external  field  '^2  are 

t)  =  i)  +  ?)},  ?)  =  'h(r,  t) 


(16) 


(17) 


and  will,  for  later  purposes,  be  regarded  as  functions  on  all  of  It  follows  directly  from  the 
definition  (5)  of  g  that 


1  C°° 

4771/-  J_^ 


1 


1 


t'  +? 


±r 
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giOdt' 
(i  -  t')^  - 


(18) 


which  depends  only  on  and  is  therefore  independent  of  the  choice  of  branch  cut. 
Furthermore,  since  (9)  ensures  that 

\t  -t'  ±f\^  \lm(t-t'  ±r)\  ^\b±q\^\b\-a>0, 


we  have 


|(f  -  t')^ -7^1  >  (1^1 -a)^  >  0  Vr. 

This  shows  that  iTi  is  real-analytic  in  K^,  at  least  if  g  has  compact  support  or  decays  sufficiently 
rapidly  to  ensure  that  the  integral  (18)  converges.  That  is,  by  taking  the  average  (17)  we  have 
managed  to  cancel  the  singularities  of  'I'(±f,  t)  on  C  as  well  as  their  jump  discontinuities 
across  V,  leaving  a  field  which  is  analytic  in  all  of  and  hence  sourceless: 

□  'I'i(f,  t)  —  0  V{r,  t)  e 


The  only  region  where  iTa  fails  to  be  analytic  is  therefore  Sa,  where  it  is  discontinuous  by 
(16).  But  even  there,  the  irregularity  is  mild  in  the  sense  that  the  jump  discontinuity 

'Ty(r,  t)  =  T'2(7,  t)  -  'l'i(r,  t)  =  t)  -  'l'(-r,  f)}  (19) 

is  bounded — unlike  that  in  any  single  branch  'Tg,  which  is  singular  on  C.  Like  'Ll  and  '^2,  'Ly 
will  be  regarded  as  a  field  on  all  of  ,  although  for  the  present  we  need  it  only  on  Sa  ■  The 
source  Ea  of  T'a^  defined  by 


477EA=n'l'A,  (20) 

is  therefore  a  distribution  supported  on  Sa-  By  the  argument  below  (10),  it  consists  of  a 
combination  of  single  and  double  layers  on  Sa ,  with  the  difference  that  these  layers  are  now 
bounded  since  Sa  avoids  the  singular  circle  C.  Even  so,  it  is  not  clear  that  the  double  layer 
can  be  realized  in  practice.  In  the  next  section  we  replace  T'a  by  a  continuous  field,  where  the 
transition  from  'Ll  to  'T2  occurs  gradually  over  a  range  of  spheroids,  whose  source  is  supported 
on  a  spheroidal  shell  instead  of  a  single  spheroid.  (This  explains  the  reason  for  viewing  'Ll ,  'I'2 
and  T'y  as  fields  on  all  of  R^.)  Such  volume  sources,  and  their  electromagnetic  counterparts 
considered  in  the  following  section,  should  be  realizable. 


2.  Spheroidal  shell  sources 


Let  H  be  the  Heaviside  step  function.  Since  0  ^  p  <  a  in  the  interior  of  Sa  and  p  >  a  m  the 
exterior,  we  have 


'L.4  =  H{a  —  +  H{p  —  a)'^2- 

It  is  natural  to  use  the  vector  field  orthogonal  to  Sp  given  by  [K3,  appendix] 

pr  +  qa 

n=yp=  2^2’ 

p‘^  + 

which  is  unnormalized  with 


2  ,  2 

2  I  ,2  P  +a- 

n  =\n\  ^  — - - 

p^  +  q^ 


^  1, 


n  •  Vq  —  0, 


V  •  n  = 


Ip 


p^+q^' 


We  can  compute  the  source  in  (20)  by  using 

V  H(p  —  a)  =  —yH(a  —  p)  —  H'(p  —  a)n  =  &{p  —  a)n. 


(21) 

(22) 


(23) 

(24) 


Applying  the  wave  operator  gives  E,4  as  a  combination  of  terms  involving  8{p  —  a),  interpreted 
as  single  layers  on  S^,,  and  S'{p  —  a),  interpreted  as  double  layers.  As  mentioned  above,  it  is 
doubtful  whether  the  double  layer  can  be  realized  in  practice,  and  this  will  get  still  worse  in  the 
electromagnetic  case,  where  the  currents  involve  one  more  derivative.  Since  we  are  interested 
in  physically  realizable  sources,  we  now  proceed  to  modify  the  above  construction.  The  terms 
involving  S{p  —  a)  and  8'{p  —  a)  are  unavoidable  as  long  as  the  source  is  conhned  to  the 
surface  Sa-  To  construct  more  realistic  sources,  we  now  choose  a  function  h  that  approximates 
the  Heaviside  function.  A  convenient  example  is 


1  1 

h{p)  =  —  arg(— p  +  ie)  =  —  ImLog(— p  +  ie),  e  >  0 

TT  7T 

which  becomes  the  Heaviside  function  in  the  limit  as  e  ^  0.  Note  that 

s 


h'(p)  = 


7t(p^  +  e^) 


(25) 

(26) 


is  indeed  an  approximation  to  the  delta-function.  But  the  simple  choice  (25)  has  the  drawback 
that  h'  is  not  compactly  supported,  so  the  resulting  source,  although  extremely  small  outside  a 
shell  of  thickness  0(s),  will  not  have  strictly  compact  spatial  support.  To  obtain  a  compactly 
supported  source,  we  now  assume  that  h  has  the  following  properties, 


h(p)  +  h{—p)  —  \  and  h{p) — 


p  ^  —s 
s 


(27) 


where  s  will  be  assumed  fixed  with 


0  <  e  <  a,  (28) 

so  that  h(p  —  a)  vanishes  on  the  disc  P  where  p  =  0.  Define  the  regularized  version  of  T'a 
by  replacing  H  by  h  in  (21), 

=  h{a  —  p)'ki(r,  t)  +  h(p  —  t).  (29) 

To  simplify  the  equations,  we  use  the  abbreviations 

fii(p)  =  fi(o' -  p),  h2ip)  ^  h(p  -  a),  (30) 

so  that 

hi^ii+h2'l>2  =  hk'l>k  (31) 


where  the  Einstein  convention  of  summing  over  repeated  indices  has  been  used.  The  source 
of  is  defined  as  usual  by 


47rE^  =□'1'^. 

To  compute  this,  note  that  since  (27)  implies  h'(—p)  —  h'(p), 
=  —h'{a  —  /7)'hin  +  h'{p  —  a)'l'2n  + 

—  h'(p  —  a)'!'  jn  + 

where  'h  j,  defined  as  in  (19)  and  given  by 


(32) 


(33) 


1 

Anir 

1 

Inif 


g{t')At' 


1  1 
Z - 1"  Z - 

t  -  t'  -  r  t  -  t'  +  r 


(t-ng(t')dt' 

(i  - 


(34) 


no  longer  represents  a  jump  discontinuity  of  the  field  since  we  are  not  confined  to  a  single 
spheroid.  By  the  same  argument  used  to  show  that  T'l  is  real-analytic  in  it  follows  that 'T j 
is  real-analytic  in  except  for  being  discontinuous  on  V  and  singular  on  C  due  to  the  factor 
1/f.  Taking  the  divergence  of  (33)  gives 

=  h"{p  —  o')'h jrP'  +  2h'(p  —  o')V'h j  •  n  +  h'(p  —  o')'hy V  •  n  + 


Since  'hy(r,  7)  is  complex-analytic  in  r  when  r  ^V, 


p  >  V'I'y  =  T'y/Vf 


where  the  prime  denotes  the  complex  derivative  with  respect  to  f , 

'I'y'  =  ^  =  Udp  +  =  ^{vl/'(f ,  t)  +  4/'(-f,  7)}.  (35) 

or  1  L 

By  (23), 

yf.n  —  {V p  —  iV q)  •  n  —  . 


Subtracting  thus  gives 

— =  h"{p  —  a)'P  jri^  +  2h'(p  —  a)'i>jfn^  +  h'{p  —  a) 'Ey  Vn  — 

But  we  have  seen  that  □T'l  vanishes  identically  and  □'I'a  is  supported  on  V,  where  /z2  =  0  by 
(27)  and  (28).  Using  (23)  gives  the  regularized  source 


»  „  T  2h'(p  —  a)  T  T 

-471 =  /;  (p-  a)'Pjn^  + - ^ +  «  )'!'/'  +  P'^ A 


supported  on  the  spheroidal  shell 


(36) 


=  {r  :  a  —  s  ^  p  ^  a  +  e).  (37) 

We  emphasize  that  E^  is  a  smooth  volume  source  that  depends  only  on  the  ‘jump  field’ 
'h y .  Taking  the  limit  e  — 0  so  that  h  becomes  the  Heaviside  function  gives  the  source  E.4 
consisting  of  single  and  double  layers  on  Sa . 


3.  Maxwell’s  equations  and  Hertz  potentials 

We  work  with  the  following  complex  combinations  of  electromagnetic  fields: 


F  ^D  +  iB  (38) 

G^E  +  \H^F-AnP  (39) 

P  =  P,  +  iP,„  (40) 

where  the  units  are  Gaussian  with  c  —  1,  Pg  is  the  electric  dipole  density,  and  P„,  is  the 
magnetic  dipole  density.  Maxwell’s  equations  take  the  form 

V-P  =  4jrp  P  +  iVxG  =  -4:Tj  (41) 

where  F  —  d,F.  In  the  general  case  of  complex  charge  and  current  densities 
P  —  Pe  tPm  J  i 

equations  (41)  are  equivalent  to 

V  •  D  —  AjtPe  i)  —  V  X  H  —  —AltJe 

V  •  B  —  Alt  Pin  P  +  V  X  E  —  —AnJm, 

so  the  imaginary  parts  (p^ ,  Jm)  represent  the  magnetic  charge-current  density.  Since  magnetic 
monopoles  are  not  observed,  we  must  require 

p,n  =  0  and  Jm  —  0.  (43) 

That  is,  equation  (41)  are  completely  equivalent  to  the  usual  Maxwell  equations  if  we  add 
the  requirement  that  (p ,  J)  is  real.  We  will  consider  solutions  derived  from  a  complex  Hertz 
potential  consisting  of  electric  and  magnetic  Hertz  vectors  [BW99,  J99] 

Z  —  Ze  +  iZfn  (44) 

whose  source  is  the  polarization, 

□  Z  =  4;rP.  (45) 

(For  this  reason,  (Z^,  Z„,)  are  sometimes  called  polarization  potentials.)  The  field  E  is  then 
given  in  terms  of  Z  by 

P  =  VxVxZ  +  iVxZ,  (46) 

and  it  follows  from  (39)  and  (45)  that 

G  =  P-nz  =  VxVxZ  +  iVxZ  +  AZ-Z 

.  ..  (47) 

=  VV  •  Z  +  iV  X  Z  —  Z. 

The  real  form  of  equations  (45)-(47)  is  [BW99,  pp  84,  85] 

□  Ze  =  47rPe  nZm  =  Alt  P,n 

E  —  VV  •  Zg  —  V  X  Z„,  —  Zg  P  =  V  X  V  X  Z;„  +  V  X  Zg 

P'  =  V  V  •  ZeS^  ^  Ze  —  Z,n  P  =  V  X  V  X  Zg  —  V  X  Z^ . 

An  inspection  of  the  expressions  for  E  and  P  reveals  the  meaning  of  the  Hertz  potentials  as 
‘superpotentials’  from  which  the  four-vector  potential  (<t>,  A)  can  be  derived  by 

o  =  -V  ■  Ze  A  —  V  X  Z,„  +  Zg. 

In  fact,  these  automatically  satisfy  the  Lorenz  condition 


(48) 


and  every  four-vector  potential  satisfying  it  can  be  derived  from  Hertz  potentials.  The  freedom 
to  choose  a  gauge  for  (<15,  A),  including  a  non-Lorenz  gauge,  is  part  of  a  much  greater  gauge 
freedom  in  (Z„  Z„)  [N55,  BW99]. 

According  to  (46),  F  is  a  curl,  so  by  (41)  the  free  charge  density  vanishes; 

p  =  0.  (49) 

Furthermore,  (47)  gives 

VxG  =  iVxVxZ— Vx.Z  =  iF, 
therefore  by  (41),  the  free  current  density  also  vanishes: 

J  =  0.  (50) 

Maxwell’s  equations  (42),  written  in  terms  of  the  microscopic  fields  {E,  B),  now  state  that 

V-B  =  0  Vx£;  +  .B  =  0 

y  •  E  —  Anph  V  y.B  —  E  —  An  Jh 

where 

-V  -Be  =  Pe  +  V  X  (52) 

represent  the  bound  charge  and  current  densities  generated  by  the  variable  polarizations 
(Pg,  Pm)-  The  fields  derived  from  Hertz  potentials  as  above  are  thus  due  entirely  to  bound 
sources^. 

4.  Spheroidal  electromagnetic  antennas 

In  this  section,  we  construct  electromagnetic  wavelets  from  scalar  wavelets  by  turning  'T  into 
Z,  then  compute  their  charge-current  densities.  It  is  essential  that  the  polarization  P  defined 
in  (45)  have  compact  spatial  support,  as  it  cannot  otherwise  be  realized.  There  are  various 
ways  to  turn  a  scalar  solution  of  the  wave  equation  into  a  vector  solution  without  increasing 
the  support  of  its  source  distribution,  the  simplest  being 

Z  =  (53) 

where  p  is  a  constant  (possibly  complex)  vector.  The  polarization  is  then  given  by 

47rP  =  pD'T  =  4:7rpE,  (54) 

so  P  and  E  have  the  same  support.  Since  E  is  a  distribution  consisting  of  single  and  double 
layers  on  V,  so  is  P.  A  similar  construction  applies  to  the  different  versions  supported  on 
the  general  branch  cut  B  and  the  spheroid  Sa-  As  explained  below  (24),  the  layers  on  B  are 
singular  on  C  while  those  on  Sa  are  bounded.  Even  so,  the  charge-current  distributions  (52) 
require  one  further  differentiation,  hence  they  generate  a  still  higher  layer  with  coefficient 
&"{p  —  a),  and  it  is  doubtful  whether  such  distributions  can  be  realized.  For  this  reason, 
we  confine  our  analysis  to  volume  sources  on  the  spheroidal  shell  (37).  Define  the  Hertz 
potential 

Z^  =  pA>l  (55) 

with  as  in  (29),  whose  polarization  density  is 

PX  =  pE^  (56) 


2 


Free  charge-current  densities  can  be  added  by  using  stream  potentials  [N55]. 


with  given  by  (36).  If  we  interpret  as  a  scalar  density,  then  (56)  suggests  an 
interpretation  of  p  as  a  (complex)  combination  of  electric  and  magnetic  dipole  moments. 
The  charge  and  current  densities  on  the  shell,  as  given  by  (52),  are 

p,  =  -Re{V.P^}  =  -Re{p.VE^}  (57) 

and 

Jfc  =  Re{P®}+Im{V  xP^}  =  Re{pE^}  -  Im{p  x  VE^}.  (58) 

Outside  the  shell  the  potential  coincides  with  Z  =  p'T,  whose  pulsed-beam  field  F 
was  computed  in  [K4]. 


5.  Extended  Huygens  sources 


F(r,  t)  = 
G{r,  t)  = 


h2{r,  t)  = 


(59) 


The  above  suggests  an  generalization  of  Huygens  sources  [Hy99],  allowing  equivalent  sources 
to  be  represented  on  shells  instead  of  surfaces  surrounding  a  bounded  source.  We  present  this 
generalization  and  compare  it  to  the  usual  method  based  on  boundary  conditions  on  an  interface 
between  electromagnetic  media.  Let  p(r,  t)  be  a  differentiable  function,  which  will  be  called 
a  zone  function.  Fix  two  numbers  pi  <  p2  and  consider  the  time-dependent  surfaces  and 
volumes  in  defined  by 

S\{t)  =  {r  :  p{r,  t)  =  pi),  82(1)  =  {r  :  p{r,  t)  =  pi) 

V'i(t)  =  {r  :  p{r,  t)  <  pi),  V2it)  =  {r  :  p{r,  t)  >  P2}. 

Given  two  electromagnetic  fields  (Pi,  Gi)  and  (P2,  G2),  with  or  without  sources,  we  want  to 
construct  an  interpolated  field  (P,  G)  so  that 

Pi(r,r),  reVlit) 

F2(r,t),  r  e  ¥2(1) 

Gi(r,  t),  r  e  yi(f) 

G2ir,t),  re  ¥2(1). 

Choose  a  differentiable  function  h2ir,  t)  such  that 

fO,  re¥x{t) 

,1,  re¥2{t) 

and  let 

h\(r,  t)  —  \  —  h2{r,  t). 

We  define  the  interpolated  field  as 

P(r,  t)  =  hk{r,  t)Fk{r,  t) 

G{r,t)  =  hk{r,t)Gk(r,t) 

where  summations  over  A:  =  1 ,  2  are  implied,  and  the  jump  field 
Fj  —  P2  —  Pi  —  Dj  +  iB  j 
Gj  =  G2-Gi=Ej+iHj. 

Then,  according  to  (41),  the  charge  density  of  (P,  G)  is 

47rp  =  V  •  P  =  Ajthkpk  +  Vfi2  •  Fj  (63) 

where 


(60) 


(61) 


(62) 


Anpk  -Fk, 


jA  =  1,2 


(64) 


are  the  charge  densities  of  the  prescribed  fields.  Thus,  in  addition  to  the  interpolated  charge 
density 

Pi  =  htpk  (65) 

we  have  a  transitional  charge  density  given  by 

47rp7’  =  Vh2’  Fj  (66) 

which  depends  only  on  the  component  of  the  jump  field  Fj  parallel  to  Vfi2-  According  to 
(60),  Pt  vanishes  outside  the  transition  shell 

V'r(0  =  {n  :  Pi  <  p(r,  t)  ^  p2}.  (67) 

Similarly,  the  current  density  is 

47r  J  =  -F  -  iV  X  G  =  AithkJk  -  hiFj  -  iVfiz  x  Gj  (68) 

where 

47TJk  —  —Fk  —  iV  X  Gk,  k  —  1,2 


are  the  current  densities  of  the  prescribed  fields.  Hence  J  is  the  sum  of  the  interpolated  current 
density 

J,  =  hkJk  (69) 

and  a  transitional  current  density  on  V7’(r)  given  by 

4jtJt  =  -h2Fj  -  iVfiz  x  Gj  (70) 


which  depends  only  on  Fj  (if  h2  is  time  dependent)  and  the  component  of  Gy  orthogonal  to 
V/!2-  The  electric  and  magnetic  transitional  sources  are  obtained  by  taking  real  and  imaginary 
parts.  Assuming  ^2  is  real,  this  gives 

4npj  =  Vh2'Dj  47Tjj  —  —k2F>j+Vh2xHj 

4Tcpj=Vh2'Bj  47Tjj——h2Bj  —  Vh2xEj. 

Letting  h2  be  complex  in  VV(f)  makes  the  transition  shell  a  chiral  medium  mixing  electric 
and  magnetic  fields.  A  further  generalization  is  obtained  by  replacing  hk  with  3x3  matrices 
(dyadics)  Mk  satisfying 


El2(r’,  t)  — 


r  e  V\{t) 
r  e  V2it)  ’ 


EIi(r,  r)  =  I  —  El2(r,  t) 


(72) 


where  I  is  the  unit  matrix.  This  makes  the  transition  shell  Vt  (t)  a  non-isotropic  medium  as 
well  as  chiral  if  EI,t  are  complex.  See  [LSTV94]  for  a  treatment  of  chiral  and  non-isotropic 
media. 

Choosing  the  zone  function  p(r,  t)  time  dependent  thus  gives  a  simple  formulation  of 
the  transition  shell  as  a  moving  source,  which  could  be  useful  in  the  analysis  of  radiation  by 
moving  objects. 

To  see  how  all  this  relates  to  Huygens’  principle,  suppose  that  we  are  only  given  a  field 
{F2,  G2)  whose  charge-current  density  (p2,  J2)  is  confined  to  Ti(f),  and  we  want  to  find  an 
equivalent  charge-current  density  confined  to  Vjit)  whose  radiated  field  in  V2it)  (but  not 
necessarily  elsewhere)  is  {F2,  G2).  We  are  free  to  choose  the  field  (Fi,  Gi)  in  any  way  that 
gives  vanishing  interpolated  sources 


Pi  —  ^kPk  —  ^I  —  ^k’^k  —  (^5 


(73) 


since  the  sources  of  the  interpolated  field  are  then  purely  transitional  and  hence  confined  to 
Vrit)  as  desired.  To  satisfy  (73),  it  suffices  to  require  that  (pi,  Ji)  be  confined  to  14(0- 


Thus,  choosing  any  field  (Fi,  Gi)  with  sources  in  14(0  and  any  function  /z2  satisfying  (60), 
an  equivalent  charge-current  density  on  Vrit)  is  given  by  (66)  and  (70).  As  pi  —>■  p2,  Vrit) 
becomes  S2  and  (pj,  Jt)  become  ordinary  Huygens  surface  sources  [Hy99]. 

The  freedom  to  choose  {Fi,  G\)  (interpreted  as  the  ‘interior  field’  if  Vi(t)  is  bounded)  and 
h2  is  constrained  by  the  requirement  that  the  magnetic  charge-current  density  must  vanish,  as 
detailed  below. 

Now  suppose  that  p  —  p(r)  is  time  independent,  so  Sk  and  I4  are  fixed,  and  choose  hk 
to  be  time  independent  and  real.  As  pi  ->  p2,  assume  that 

lim  V/!2(r)  =  S{p(r)  —  p2)n(r) 


(74) 


where  n{r)  is  a  normal  vector  field  on  ^2  pointing  into  14.  Then  (66)  and  (70)  give 
PT  &(p{r)  -  P2)(y 

Jt  &{p(r)  -  P2)K, 

where 

Ana  —  n-Fj  and  An K  — —in  x  G j  (75) 

are  the  surface  charge  and  current  densities  on  S2,  whose  real  and  imaginary  parts  give  the 
electric  and  magnetic  surface  sources: 

47rCT^  =  n-Dj  AnK^  =  nx  Hj 
Ana,„  —  n  -  Bj  AnK,n  =  — n  x  Ej. 

Since  magnetic  monopoles  are  not  observed,  and  must  vanish.  When  pi  <  p2,  this 
may  be  accomplished  if  h2  can  be  chosen  so  that 

V/Z2  •  By  =  0,  Vh2  X  By  =  0,  (77) 

which  is  possible^  if 

By.By=0  WreVrit).  (78) 

In  the  limit  pi  —y  p2,  (77)  reduces  (76)  to 

n-Dj  =  Anae  n  x  Hj  =  AnKe 

n  •  By  =  0  n  X  By  =  0,  ^  ^ 

which  are  the  usual  boundary  conditions  on  an  interface  between  two  media. 

Returning  to  the  general  time-dependent  setting,  consider  now  an  alternative  procedure  of 
special  interest  here.  Instead  of  interpolating  two  prescribed  fields  (Fh,  Gt),  let  us  interpolate 
two  Hertz  potentials  Zj : 

Z^hkZk.  (80) 

As  seen,  this  automatically  results  in  vanishing  ‘free’  sources  p  =  0,  J  =  0.  The  polarization 
is  found  to  be 

471 P  —  \Z\Z  —  An hj^Pji  -t  2A2Z y  —  2(V/22'V)Zy  —  ZyA/72 
where  AnPt  =  OZh  are  the  polarizations  of  the  prescribed  fields.  If 
suppBi  C  14  and  suppB2  C  Vi, 

then  the  interpolated  polarization  hkPk  vanishes  and  the  polarization  is  purely  transitional  on 
Vrit): 

AnP  =  2k2Zj  -2(Vh2-y)Zj  -  ZjAh2.  (81) 

This  generalizes  (56),  and  the  bound  charge-current  densities  derived  from  P^  and  B,„  via 
(52)  generalize  (57)  and  (58). 

^  Letting  Vt  be  a  non-isotropic  medium  by  using  Hyt  (72)  makes  it  easier  to  enforce  the  absence  of  magnetic 
monopoles. 


6.  Conclusions 


We  have  improved  on  the  computation  of  sources  for  electromagnetic  wavelets  given  in  [K4] 
in  two  ways:  (a)  the  spheroidal  surface  Sa  supporting  the  sources  has  been  replaced  hy  a 
spheroidal  shell  supporting  smooth  volume  sources.  This  eliminates  the  multiple  layers  on 
Sa  which  make  the  sources  difficult  if  not  impossible  to  realize;  (b)  by  deriving  the  sources  from 
Hertz  potentials,  we  have  eliminated  the  magnetic  charge-current  density,  further  facilitating 
their  realizability. 

The  problem  with  the  magnetic  sources  in  [K4]  can  be  better  understood  from  the  current 
perspective.  Let  the  zone  function  be  p  =  Re  f  and  Sk,  Vjt  be  as  above  with  0  <  pi  <  p2-  Let 
the  fields  (Fk,  Gk)  be  derived  from  the  Hertz  potentials 

Zk  =  P'^k,  k  =  1,  2 

with  'I'i:  given  by  (17).  Recall  that  'I'k  are  analytic  in  (f ,  t)  for  p  >  0.  Therefore  the  fields 
{Fk,  Gk)  are  analytic  in  (r,  t)  for  r  ^V,  and 

Pk  =  0,  Gk^Fk^  Ek+iBk  Vr  ^  V. 

The  jump  fields 


Fj  —  Gj  —  E]  +  \Bj 

are  also  analytic,  as  is  their  polarization  scalar  [K3a] 

F]  =  Fj.Fj^E]-B]+2\Ej.Bj.  (82) 

The  condition  (78)  thus  requires  the  imaginary  part  of  an  analytic  function  to  vanish  for  r  e  Vt, 
which  implies  that  Ej  vanishes  identically.  Although  every  electromagnetic  field  must  have 
Ej—kOin  the  far  zone  [B15],  fields  satisfying  Ej  =  0  globally,  called  null  fields,  are  rather 
degenerate.  In  particular,  the  electromagnetic  wavelet  fields  are  not  null  and  hence  cannot 
fulfil  (78).  Instead,  we  have  begun  with  an  interpolated  Hertz  potential  (80)  and  derived  bound 
sources  in  the  transition  shell,  thus  preserving  analyticity  without  invoking  the  existence  of 
magnetic  monopoles. 

There  is  still  an  unsatisfactory  aspect  to  the  polarization  (56)  and  charge-current  density 
(57),  (58).  Namely,  they  depend  on  the  fixed  vector  p  and  thus  do  not  conform  to  the  spheroidal 
geometry.  This  suggests  using  methods  of  constructing  Z  from  'T  other  than  (53).  While 
(53)  is  the  complex  version  of  Whittaker’s  potentials  [W4],  there  are  alternatives  which  do  not 
require  a  fixed  polarization  vector,  such  as  Debye  potentials;  see  [BD98]  for  example.  Such 
alternatives  will  be  considered  in  future  work. 


Acknowledgments 

It  is  a  pleasure  to  thank  Drs  Richard  Albanese,  Grant  Erdmann,  Sherwood  Samn  and  Arthur 
Yaghjian  for  helpful  discussions.  I  am  also  grateful  to  Dr  Arje  Nachman  for  his  sustained 
support  of  my  research,  most  recently  through  AFOSR  grant  no  FA9550-04-1-0139. 

Note  added.  After  this  paper  was  finished,  I  learned  from  Dr  Arthur  Yaghjian  that  a  similar  formulation  of  the  gradual 
transition  between  two  electromagnetic  media  has  been  developed  by  Lindell  et  al  [LTNO].  The  present  form  is 
somewhat  more  general  in  the  following  respects,  (a)  The  regions  and  surfaces  are  allowed  to  be  time  dependent, 
(b)  The  derivation  of  generalized  Huygens  sources  in  [LTNO]  makes  the  assumption  that  one  of  the  prescribed  fields 
vanishes,  which  we  find  to  be  unnecessary.  This  allows  extended  Huygens  sources  on  Vj  with  nonvanishing  ‘internal’ 
fields. 
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